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ON ' We investigate the possible origin of anomalous Hall effect in the CMR (colossal magnetoresis- 

tance) materials - the doped rare earth manganites - observed recently by Matl et al It is 
'_ , demonstrated that the spin-orbit interaction in the double exchange model couples magnetization 

to the Berry phase associated with three dimensional spin textures and induces a non-zero average 
topological flux which in turn generates an anomalous contribution to transverse resistivity. The 
same effect, but involving the orbital Berry phase, occurs in the model with orbital degeneracy and 
Coulomb repulsion. 
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It has been known for a long time ^ that the Hall effect in ferromagnets exhibits a low field "anomaly" which can 
be parametrized as magnetization, M , dependent contribution to transverse resistivity: p xy = RhB + RaM. The 
anomalous Hall constant Ra is usually explained in terms of skew scattering due to spin-orbit coupling as reviewed in 
Ref J2|. Recently Matl et al Q| measured the Hall resistivity of Lai-^Ca^MnOs with x = 0.3 and found negative Ra (in 
contrast with positive, hole like, Rh) which becomes significant above \QQK and peaks about 20-ftT above the Curie 
temperature T c w 265K. This temperature dependence (which follows closely that of the longitudinal resistivity) 
| contrasts with the behavior of, for example, Ni, where Ra peaks well below T c and contradicts to the existing theories 
P,pp5|]. There exist also related Hall measurements on similar material at relatively high temperatures Q and on 
La^-j-Ca^CoOs (0.1 < x < 0.5) 0. These new data lead us to revisit the theory of the anomalous Hall effect. 

The Manganites are described by a model which includes double exchange, Jahn- Teller, and electron-electron 
(Hubbard) interactions: 



^ I Hq = —t ^ cl v gCr+^a - J H ^ S r • c\, ^ a T acr 'Cr,(T',a + #JT + #e-e ■ (1) 
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where aJ CTa is the creation operator for the conduction electron with spin a and the orbital index a = 1,2, t ~ 
0.1 — 0.5eV is the hopping matrix element and S r is the core spin. This model is suggested by the electronic structure 
55 ' of the doped rare earth manganites |]8|-|l0| (Rei-^A^MnOa where Re is a rare-earth element such as La or Nd and 
0^ . A is a divalent metal ion such as Sr or Ca) where 1 — x electrons occupy the orbitally degenerate e g states (d x 2_ y 2 
' and d 3z 2_ r 2) and couple strongly to the core S = 3/2 moment via Hund's rule, Jh ~ 1 — 2eV. The on-site Coulomb 
repulsion B e - e = UJ2 r n r,i;,a- n r,l,a + VJ^r.aM 71 ^, 1 ^^ (where n r ^^ a = cl a a c r ^^) is believed to be of the same 
order as Jh in addition, it has been argued [ p^|Jl4| that in the rare earth manganites an important role is 

^ , played by the Jahn- Teller (JT) coupling of the e g electron with the phonon mode Q describing a local tetragonal 
distortion which lifts orbital degeneracy via Hjt — g^2 r c\. lj a Qlf'c r . at b- This coupling can explain the transition from 
metallic to activated transport at T c for a certain range of doping x. Below we shall discuss both the case of large 
O ■ local distortions where the orbital degree of freedom is essentially quenched and the case where distortions are weak 
, and the orbital degeneracy is retained (at least on the scale of ksT) which may be appropriate for the ferromagnetic 
metal phase of the compounds in question. 

Since the anomalous Hall effect manifests itself as the dependence of p xy directly on the magnetization, it must 
originate from the spin-orbit interaction Q. It includes the (relativistic) interaction of the spin of the itenerant 
electron with the local electric field and the magnetic interaction of the itinerant electron with the core spins. The 
latter term, S • p x r/r 3 , because spin-orbit coupling is quenched by the cubic crystal field has no matrix element 
bewteen e g electrons on the same site, but it does have matrix elements bcwtccn different sites. In the continuum 
limit, the sum of both contributions can be written as lilMffl 



H so =i J d r u ab e ijk Si(r) djcl a (r) d k c a ^{r) . (2) 

As with the hopping matrix, we shall for simplicity ignore the orbital dependence of the coupling constant and 
take u a b = S a bU. With the spatially uniform spin-orbit interaction mjn the coupling to the uniform component of 
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S(r), i.e., magnetization, has the form of a total derivative and can be integrated out to the boundary. Yet, we 
shall now show that in the presence of interactions this topological term is important and generates the anomalous 
Hall effect. Previous authors have studied spatially non-uniform spin configurations due to magnetic impurities or 
spin-wave excitations about the ferromagnetic ground state. Here we argue that Berry phase effects caused by non- 
coplanar spin configurations give the dominant contribution to p xy , at least in the limit of strong carrier-spin couplings 
appropriate to the maganites. 

As a mathematical convenience to make the Berry phase physics explicit, we use the standard Hubbard-Stratonovich 
approach in which we express the interacting electron in terms of a non-degenerate spinless fermion ip r and the slave 
boson fields z Tt(J and w rya which parametrize the spin and the orbital degrees of freedom respectively: c CT , a = xj^ z a w a - 
The bosons obey the single occupancy constraint: z r ,oZ r ,<j = ^2 a Wr,a w r,a = 1 an d the average fermion density is 
(■i/^Vv) = x. Neglecting terms of relative order £/Jh and t/U, we find 

r,[i r r 

where = (£f , z^j and = (wi, Ub). The factorization of c CTja not only enforces single occupancy but projects onto 
the low energy spin and orbital state. Assuming Jh > t the carrier spin and hence s = zjrz r is aligned with S r . 
The generally complex spin overlap factor z£_|_„z r is then determined by the Euler angles of the corresponding local 
moments and is identical to the one found in a conventional derivation of double exchange fiof . Orbital degree of 
freedom results in an analogous orbital overlap factor wj_|_„w r . Note however that since the tetragonal/orthorhombic 
phonon field Q = Q\T X + iQ2T v is real, large g\Q\ forces the orbital "isospin" n = w^fw into planar configurations. 
We emphasize that although this formalism is the most natural one in the strong coupling limit, we expect the physics 
to apply at any coupling, with possibly different amplitude as discussed below. 

The Hamiltonian (|^) is invariant under two local £7(1) transformations JlJ,|l^]. One of them involves the spin degree 
of freedom and is given by ijj r — > ip r e l ^ r and z r — * z r e l ^ r . This spin i7(l) symmetry is broken in the ferromagnetic 
(FM) phase, but is restored in the paramagnetic (PM) phase. The other one involves the orbital degree of freedom 
and represented as tp r — > i/ve i( = r and w r — > w r e'^. As long as the orbital ordering does not occur Jl3|,[l2| the orbital 
U(l) symmetry is preserved. These U(l) symmetries give rise to two gauge fields which are defined as 

a t = ^[(^^z-zt^z] 

<=|[(^wt)w-wt^w] . (4) 

Non-trivial configurations of the gauge field are generated by non-coplanar arrangements of spin or orbital moments; 
they affect carrier motion just as does an external magnetic field, and may therefore give rise to an anomalous Hall 
effect. Their coupling to the charge, spin, and orbital degrees of freedom is particularly transparent in the continuum 
limit HjHJ of (|): 

tfcont = J d 3 r [ieff^^ + eA„ + a* + a£) 2 V> + ^(ify + a* ) 2 z + pSffW^ + a^) 2 w + g^Q(r)w] (5) 

where t e s is the effective hopping parameter of charge carriers and the spin and orbital stiffness are defined in the 
sense of a self-consistent mean field theory. Q(r) is the quenched random phonon field. The scale of p s c 'g is set by the 
electron kinetic energy while the effective hopping of the holes at finite temperatures is reduced from its bare value t 
by the fluctuation of the local moments and the orbital degree of freedom. The (tpip') = 1 — x factor appearing in the 
MFT, can be absorbed into resetting the boson density z^z = w^w = 1 — x. We have introduced the electromagnetic 
gauge field A which only couples to the charge of the fermion (e > since ^ creates holes). The key point in (|5|) 
is that the fermions feel an effective magnetic field B e g = B + e~ 1 b s + e~ 1 b™ = V x (A + e _1 a s + e _1 a") which 
includes the contribution of the bosonic Berry phases described by & s,n . 

The spin-orbit Hamiltonian (|^) includes the coupling to the Berry phase which is made explicit by substituting 

Ca,a = Ip^Zo-Wa. 

H^ D rj J dr 3 u{^) M • V x (a™ + a s ) (6) 

where we have assumed that z and w are slowly varying functions and suppressed the terms with derivatives acting 
on ip. We observe here that the spin-orbit interaction couples the average magnetization to the bosonic Berry phases 
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(topological flux) b s '". The latter arises through non-planar configurations of the s(r) and/or n(r) via bf (r) = 
CyfeS • djs x <9fcS and analogously for b". Note however that strong Jahn- Teller effect (i.e. g |Q| 3> fcflT) would confine 
n to a plane and suppress b™. The fact that the spin-orbit interaction couples magnetization to the topological flux 
associated with spin (or orbital) textures is quite general and not restricted to the strong coupling limit and any slave 
boson/fermion parametrization. However, the effects become weak for small Ju/t. 

Let us start by considering the case of quenched orbital phases, b™ — and analyse the effect of spin textures in 
the paramagnetic and ferromagnetic phases. In the FM phase z bosons are Bose-condensed corresponding to FM 
order in s. As a result, the spin gauge field a s becomes massive with the mass proportional to the density of the 
condensate and thus there is no uniform contribution to b s . On the other hand in the PM phase we expect to find a 
quadratic contribution, (b s ) 2 , to the total energy and a finite topological susceptibility defined as \t = d(b s ) / dM . 
The calculation proceeds by integrating out the ip and z fields in Eq.([5]) and results in an effective Hamiltonian: 

H b = J d 3 r [ XF (eB + b s ) 2 + X b (b s ) 2 + u(l - x) M • b s ] (7) 

with the diamagnetic susceptibilities xf = K~ 2 t 2 s D(eF) and xb — n (0)y / 2Tp^ s (where n(e) = 1/ (e^^ B ^ T — 1) is 
the occupation number of the z bosons) are determined by the current /current correlators of the fermions and bosons 
respectively. Minimizing with respect to b s yields 

(b')=- J (1 - X) M-^—B (8) 
2{ X f + Xb) Xf + Xb 

from which we read off the topological susceptibility xt = — 2(xf+xb) ' Magnetization here refers to the average local 
moment in units of g e p,BS where g e is the Lande g-factor and pis is the Bohr magneton. The topological flux per 
plaquette of the lattice is of order u/xf.b which we can estimate conservatively by taking u ~ 1 — 10 A" and x ~ * 
yielding (b s ) ~ 10~ 3 — 10~ 2 Z~ 2 where I is the lattice constant. This flux is comparable to eB/hc for IT field. One also 
observes that the flux is maximized for a nearly full band as may be expected for the effect arising from interactions. 

According to the Ioffe-Larkin phenomenology [pl| the constraint between the fermionic and the bosonic currents 
which follows from (|J) implies that the Hall resistivity p xy is given by the sum of fermionic p xy and the bosonic p xy 
contributions. Using the semiclassical Hall constants corresponding to fermion density np — x and boson density 
tib = 1 — x, we have p xy — (eB + b s )/e 2 np and p xy = b s /e 2 ns so that 

u(l -x)(n f 1 + n B 1 ) XBnp 1 - XftIb 1 B 

Pxy = 7TT7 i \ M i • ( 9 ) 



2e 2 (xF + XB) Xf + Xb 



from which we identify: 



Rh 



1 Xsrip 1 - XFn B X 

e Xf + Xb 
u(l - x^rip 1 + n B ') 



Ra = 9 2 * \ ( 10 ) 

2e 2 (xF+XB) 

According to the naive classical estimate, u > and Ra has the sign ^| opposite to Rh- This could explain the 
sign of the anomalous Hall effect observed by Matl et al M. Note also the compensation of Rh by the bosonic 
contribution: the large apparent carrier concentration observed in the experiment suggests that some compensation 
effect does indeed exist. However, the calculation of the Hall constant in an interacting system clearly requires greater 
effort than the semiclassical estimates. 

The above estimate for xf given above is appropriate for the metallic paramagnetic phase (such as observed in the 
doped rare earth manganites at the values of x somewhat higher than optimal value for the CMR effect); at lower 
dopings the PM phase has very high resistivity, thus we expect in this situation that xf "C xb [D(€f) — * 0) and the 
relevant carrier density np ( rip 1 ^s> rig 1 ) to be thermally activated. For the CMR compounds, the insulating PM 
phase is likely to be caused by the large local Jahn- Teller distortions (l^Jl^l which split the conduction band into 
filled orbitally non-degenerate subband with spectral weight 1 — x and an empty orbitally degenerate subband with 
spectral weight x. Conduction is due to the thermally excited carriers in both subbands and the proper determination 
of the sign and magnitude of Rh is non-trivial. (Alternatively one may think of transport in terms of small polarons 
||.) However the appearance, due to spin-orbit interaction, of the average topological field (b s ) as well as its coupling 
to the carriers are robust. We expect Ra ~ e~ 1 xrRH with xt = —u(l — x)(2xb)~ 1 , which follows from (10) under 
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Xf ^ Xb and rip rig 1 assumptions but can be understood simply by saying that (b s ) acts just like magnetic 
field. 

Let us now discuss the ferromagnetic phase. As already noted, in the FM phase the a s gauge field acquires a mass 
due to the divergence of \b at fc = 0. This phenomenon is reminiscent of the Meissner effect in a superconductor, 
but there is a crucual difference due to the fact that the topological flux on the lattice is not conserved. A continuum 
FM spin texture with b|(r) — e^feS • djS x dkS is the Belavin-Polyakov |23| skyrmion which in a 3D ferromagnet 
would be a line defect (akin to a flux tube) with the line tension 2-7T/9g ff |q| where q is the topological charge equal to 
the integral of b s over the crossection of the "tube" divided by 87r. Yet when this continuum solution is transfered 
to the lattice, the quantization of \q\ is lost since shrinking the lateral size of the skyrmion allows the texture to "fall 
through" the plaquette. For the same reason (and in contrast to a flux line in a superconductor) the skyrmion line can 
have free ends and at T ^ there will be a finite density of such objects. Alternatively one can argue that the spins 
around any plaquette may be deformed into a texture carrying topological flux q — at an energy cost a/9g ff |g| (a 
being a constant of O(l) ) which should be contrasted with the quadratic energy cost in the PM. One can calculate 
the topological susceptibility from the partition function Z{M) = J dh s exp[— (ap^ ff /87r)|b s |/r — uh s ■ M]. One finds 
Xt ~ —8iTuT/ap^ s for T / p s cS <C 1. Since the topological b couples to the carriers exactly like the magnetic field we 
expect: Ra — 6~ 1 \tRh- Note that approaching T c the topological susceptibility \t should be enhanced both by the 
softening of p s eS (due to the reduction of electron kinetic energy Jl3| which underlies the FM exchange constant) and 
by the contribution of the large scale skyrmion textures which take advantage of the vanishing spin stiffness at long 
wavelength. The enhancement of \t and hence Ra could explain the observed upturn of Ra below T c . 

While the temperature dependence of Ra indicated by the above arguments is qualitatively consistent with that 
observed by Matl et al 0] the region close to T c requires a more careful treatment. At present we neither understand 
why (or whether) Ra(T) ~ p(T) nor why both quantities peak 20 — 30if above T c . An a priory calculation of Rh is 
beyond our present scope. However we did demonstrate how in the presence of interactions the spin-orbit coupling can 
generate direct dependence of p xy on magnetization via the appearance of topological flux. This should be contrasted 
to the traditional view MM] of skew-scattering contributing to p xy in the 3d order via ((M - (M)) 3 ) which peaks well 
below T c . 

Curiously, the analysis presented for the PM phase applies equally well to the description of the orbital Berry phase 
effects in the orbital "liquid" state [jl2| (.g|Q| <C fcsT) where the orbital degeneracy is restored in FM phase. While 
in double exchange problem the Hubbard U played no role as the single occupancy was ensured by the Hund's rule 
splitting, if the orbital degeneracy is restored, the effect of U becomes important and necessary for the appearence 
of the topological effects. In this case, b s in Eq.^ and Eq.|| is replaced by b™ and the rest of the discussion goes in 
the same way. Therefore, if U is large and g|Q| small, the effect of the orbital Berry phase would provide additional 
contribution to the anomalous Hall coefficient in FM phase. 

In conclusion we emphasize the generality of the physical mechanism by which the topological textures in magnetic 
systems get polarized by the magnetization via spin-orbit coupling and contribute to p xy . Comparing an itenerant 
magnet to the double exchange magnet considered above, the former corresponds to a weak coupling regime, where 
the local moments and hence any possible textures with (S; • Sj x S^) ^ disappear at T c . By contrast, manganites 
with their core spins and large Jh are in the strong coupling limit. While our work points to the possible origin of the 
unusual anomalous Hall effect in doped manganites much further work is required to make a quantitative comparison. 
It would also be interesting to measure the anomalous Hall effect for Lao.7Sro.3Mn04 or the like compound, which 
remains metallic in the PM phase, which could be compared with our result on PM metallic phase. 
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